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1. Introduction 



The LHC will soon resolve the question as to whether low energy supersymmetry is 
the solution to the hierarchy problem; and if it is, moreover, the LHC and a future e + e~ 
linear collider (LC) will lead to very precise measurements of the sparticle spectrum and 
couplings. The success of gauge unification in the MSSM suggests a Desert, the existence 
of which would mean that extrapolation of the MSSM couplings and masses to high scales 
will lead to immediate information about the underlying theory; for example regarding 
the commonly assumed universality of soft scalar masses, gaugino masses and cubic scalar 
interactions. 

One component of this analysis is the running of masses and couplings between the 
weak and gauge unification scales, which is governed by the renormalisation group (3- 
functions. In this paper we compare the results for this process using one, two and three- 
loop /3-functions. In each case we generally use the same one-loop corrections for the 
relationship between running and pole masses for the various particles, with some use of 
two-loop results such as for the top quark mass. We anticipate that by the time sparticles 
are discovered complete two-loop threshold corrections will be available; the effect of these 
we would expect to be of the same order of magnitude as the effect of using the three- 
loop (as opposed to two-loop) /3-functions, which, as we shall see, is surprisingly large for 
squarks. 

The plan of this paper is as follows. In section 2 we review the exact results that 
relate the (3- functions for the soft masses and interactions [|1] -0 to the /9-functions of the 
dimensionless gauge and Yukawa couplings Q- 0, which we then give through three loops 
for the MSSM generalised to incorporate and nio sets of SU$ 5(5) and 10(10) represen- 
tations respectively. (A motive for grouping additional matter in this way is that complete 
SU5 representations do not (at one loop) change the prediction of sin 2 9w (or alternatively 
of g%(Mz)) that follows from imposing #1,2,3 gauge unification. Also unchanged at one loop 
is the gauge unification scale, Mj; but at higher loops this scale increases and can ap- 
proach the string scale.) We also give a simplified example of a three- loop soft (3- function; 
general results for all the f3- functions are available at Ref. J7|. 

In section 3 we present and discuss our results for the sparticle spectrum for a set 
of Snowmass Benchmark Points 0], all corresponding to the standard universal boundary 
conditions at unification, except for one case with non-universal gaugino masses. We 



compare our results with the useful website Ref. B (see also Refs. [pjj , |Tl| 



In section 4 we consider the effect of additional matter fields in SU$ representations, 
as discussed in Refs. |TjJ (for earlier work see for example Refs. jrjj) and by ourselves 
in a previous paper [ fL5| . We give some further examples of the effect on the sparticle 
spectrum of such matter. Finally section 5 contains our conclusions. 



2. The Soft Beta functions 

For a general N = 1 supersymmetric gauge theory with superpotential 

the standard soft supersymmetry-breaking scalar terms £1X6 cLS follows 
V soit = (! 6<i && + \h ijk ^ 3 4> k + c.c.) + (m 2 Y 



(2.1) 



(2.2) 



where we denote 4> l = 4>* etc. (For the generalisation to the case when Vg ^ includes a 
term linear in (p see [|16| . ) 

The complete exact results for the soft /3-f unctions are given by: 
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where 7 is the matter multiplet anomalous dimension, and 
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Here M is the gaugino mass and Y ijk = (m 2 yiY jkl + (m 2 ) J \Y ikl + (m 2 ) k iY ijl . Eq. ( pO| ) 
holds in a class of renormalisation schemes that includes the DRED'-one[|l7[], which we will 
use throughout. 

Finally the X function above is given (in the NSVZ scheme |18[]) by 
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16tt 2 [1 -2g 2 C(G)(167T 2 )- 1 ] 
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(2.5) 



where 

S = r-Hr[m 2 C(R)} - MM*C(G), (2.6) 

C(R),C(G) being the quadratic Casimirs for the matter and adjoint representations re- 
spectively. There is no corresponding exact form for X in the DRED' scheme []T7[|; we will 



require the leading and sub-leading contributions, which are given by ||19|| : 

167r 2 X DRED '« = -2g 3 S (2.7) 

and 

(167r 2 ) 2 X DRED '( 2 ) = (2r)- 1 g 3 tv[WC(R)] - 4g 5 C(G)S - 2g 5 C(G)QMM*, (2.8) 

where 

W j i = )-Y ipq Y™ n {m 2 y n + \Y^Y pqn {m 2 ) n i + 2Y ipq Y^ \m 2 )\ 

£ ^ (2-9) 

+ h ipq h jpq - 8g 2 MM*C(R) J i. 

and Q = T(R) — 3C(G), and rT(R) = tr [C(R)], r being the number of group generators. 

We now present the results for the gauge /3-functions and anomalous dimensions. 
These results are valid in the DRED' scheme |jT7| (or indeed the DRED one|2(J, which 
differs from DRED' only when we come to the soft /^-functions) . The MSSM superpotential 
is: 

W = H 2 QY t t c + HiQYtb + H\LY T r c + iiH x H 2 (2.10) 
where Y t , Yt,, Y T are n g x n g Yukawa matrices i , and we define 

T = Y t Y}, B = Y h Y h \ E = Y T Y},f= Y?Y t , B = Y^Y b , E = Y?Y T . (2.11) 

The SU3 ® SU2 ® U\ gauge /3-functions are as follows: 

(3 9i = (IGtt 2 )- 1 ^ 3 + (16tt 2 )- 2 ^ 3 M£M " + ( 167r2 )" 3 ^ } + • • • ( 2 - 12 ) 

where 

b% = hn 5 + |n 10 + h = hn 5 + §m + 1, 63 = hn s + %n 1Q - 3 



ai = ^trT + ^trS + ™tiE, a 2 = 6trT + 6trS + 2tvE, a 3 = 4trT + 4trS 



(2.13) 



Vt,b, T here are the transposes of the Yukawa matrices used in Ref. ||. 
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and 




(2.14) 



For the anomalous dimensions of the chiral superfields we have at one loop: 



(167r 2 ) 2 7 f } = -2T 2 - 6(trT)T - 2Y? BY t + (6a 2 - fa 2 ) f 

+ (ft&i + + Wdh! + (lbs + f)9l (2.16a) 
(167r 2 ) 2 7 f } = -2B 2 - 6(trS)S - 2Y$TY b - 2(trE)B + (6g% + f # 2 ) B 

+ (¥i + 4-M + §9hi + (lbs + fhi (2.166) 

(167r 2 ) 2 7 g } = -2T 2 - 3(trT)T - 2S 2 - 3(trS)S 

- (trE)B + a 2 (f T + \B) + ±g 2 l9 2 + 8g 2 g 2 + ^g 2 g 2 

+ (lbs + f)9 4 s + (f 6a + f )a 2 4 + (±h + ^)gi (2.16c) 

(167r 2 ) 2 7 g 2 } = -9trT 2 - 3trST + (l6o 2 + fa 2 ) trT + (f b 2 + § )g\ 

+ ^2 2 + (^i + T§o)^ ( 2 - 16rf ) 
(167r 2 ) 2 7 g i ) = -9trS 2 - 3trST - 3(tr£ 2 ) + (I603 ~ ff 1) trS + hl tvE 

+ (§&a + !)^2 4 + fefol + + ifo)<rf, (2-16e) 
(167r 2 ) 2 7 i 2) = -2E 2 - E(3trB + trE - fa 2 ) + (f b 2 + |)a 2 4 

+ ^2 2 + (^i + i§o)^ ( 2 - 16 /) 

(16tt 2 ) 2 7 | 2 ) = -2E 2 - E(6trB + 2tvE - 6g 2 + fa 2 ) + (f h + § )g\. (2.16a) 
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(2.15) 



and at two loops [pi 



The three-loop gauge /3-function terms are (henceforth we suppress (167T 2 ) L factors). 



p(f =g f \ftrT 2 + 18(trT) 2 + f trB 2 + f (trB) 2 + f trTB + f trT 2 + f (trT) 2 
+ ftrTtrT - (if g\ + fg 2 + *§g 2 ) trT - (f g 2 + fg 2 + 2 -§g 2 ) trB 
- (M + fti^E + (±f - In 2 - ( ™ + 6n 10 )n 5 - l§Vo - fn 2 10 )gl 

( 27 2 I / 27 i 9„ \„ , 81 , 261 „ , 27„2 \„4 / 1 , 27„ , 69\ J„2 



( 7 J2 , / 7507 , 9_, \„ , 12859 ^ 
liO n 5 + ("900" + 4 n 10j«5 + -300~ n 10 



207 2 , 32117 \ 4 
40 ' J 10 + 375 >9l 



(2.17a) 



=£ ,| 24(trT 2 + trT 2 ) + 18[(trT) 2 + (trB) 2 ] + 12trTT + 12trTtrT + 8trT 2 

+ 2(trT) 2 - (32s 2 + 33s 2 )(trT + trB) - g 2 (f trT + ^trT) - (llg 2 + fg 2 ) trT 
- [(6ni + 18)n 5 + ^Vo + 18n 2 - 4%| + [(8n 10 + 24)s 2 2 - (^n 10 + |)^ 2 ]^ 
" (f™5 + (f + f ™io)n 5 + f nio + ^n 2 - 35)^ + (±n 10 + § + ^n 5 )<7 2 <7 2 2 



f9„2 | /3_, , 441\_, , 9 „2 i 457 , 1513 „ \_4~ 
~ Ii0 n 5 + U n 10 + 100 J n 5 + 40 n 10 + ^5" + "300~ ni °J#lJ ' 

/?g) =#3 18(trT) 2 + 12trT 2 + 8trTT + 18(trT) 2 + 12trT 2 + 6trTtrT 

- C-f9 2 s + l2g 2 )(trT + trB) - 2 (ff trT + fftrT) 

+ (MI + 2 -fn 10 - ^n§ + ( 2 f- fn 10 )n 5 - fn 2 )gi 
+ [(2ni + Q)gl + (§nio + 3^5 + D^ 2 

- [(¥ + l n io)n 5 + f n 2 + ^n 10 + 27]^ - (|m + f )^| 

/ 1 „2 , / 2689 I 3„ \„ i 27„2 , 1702 , 3353 „ \-4 
~ Kw n 5 + 1"900 + 4 n W) n 5 + 20^10 + ^f5~ + "300" ni 0^1 



(2.176) 



(2.17c) 



The three-loop results for the anomalous dimensions are as follows: 
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7 g) = k(T 3 + T 3 ) + ABTB + ATBT + 6T 2 trT + T 2 (6trT + 2trT) 

+ B[6tv(TB) - 9(trT) 2 - 6trTtrT + 18tr(T 2 ) - (trT) 2 + 6tr(T 2 )] 

- 9T(trT) 2 + 18Ttr(T 2 ) + 6Ttr(TT) + g 2 [T 2 - k) + B 2 - ±fc)] 

+ [(3/c - 3)^ 2 + fg 2 3 ](T 2 + B 2 ) + [(2 - f % 2 + 18^ 2 + (8k - 8)g 2 ]TtrT 

+ gfB[(f - §/c)trT + (§/c - §)trT] + # 2 T(18trT + 6trT) 

+ 8<7 3 2 T[(/c - l)trT + trT] + gtT(^k - - §n 10 - gn 5 ) 

+ ^dfe* - ff - gmo - &TI5) - (|§trT + ^trT + ^trT),? 4 

+ (§* - D^lT + - %)g\g\B + g 2 g 2 T(§k - f ) + ^(g* " D 

- ^ 4 [(T + T)(f + f k + f n 10 + |n 5 ) + f (trT + trT) + f trE] 

- 4g 2 g 2 (T + B) + gj[(T + T)(§ - if k - 12n w - 4n 5 ) - f (trT + trT)] 
+ (f " 35*^10 - ^kn 5 ~%k + + ^n b )g\g\ - \g\g\g\ 

+ 0i(?ls55 ~ fc (^™io + T^o n 5 + S) + (^o n 5 + § + ^10)^10 + i^™5 + iio n 5) 

+ - ^ fallO - 255*^5 - + A^IO + 4^5) 

+ ^ 4 (ff - |fcn 10 - £fcn 5 - ff k + fn 10 + fn 5 ) 

+ dtdii^ ~ 6kn 10 - 18k + 24n 10 - 2n 5 ) + glgj(8 - 4kn 10 - 12k + 14n 10 - 2n 5 ) 
+ gU^f + fkn 10 + fkn 5 + if k + |n 10 n 5 + f n 10 + fn 2 l0 + f n 5 + |n 2 ) 
+ 9l( 2 ~W + Hfnio + fn 5 + if ) + n 10 (4n 5 + 2 -f + 6n 10 ) + f>n 5 + |n 2 ) 

(2.18) 



7 [ 3) = fci? 3 + T 2 (6trT + 2trT) 

+ T[6tr(TT) - 9(trT) 2 - 6trTtrT + 18tr(T 2 ) - (trT) 2 + 6tr(T 2 )] 

+ # 2 T 2 (9 - §fc) + (8 - 2% 2 TtrT + (3/c - 3)#fT 2 + # 2 T(18trT + 6trT) 

+ (8k - 32)g 2 EtrB + g\E(-^§ + ^k - f n 10 - §n 5 ) (2.19) 

+ <7 4 (-ftrT-fitrS-fJtrT) 

+ 9l9lE(~% + %k) + gtE(-f -fk- 2 in 10 - f n 5 ) 
+ <7!(-f trT - f trT - f trT) + S 
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7 < 3) = (6 + 2k)f 3 + 6T 2 trT - 2Y} BTY t - 2Y?TBY t + 6Yj B 2 Y t 

+ T(36tr(T 2 ) + 12tr(TS) - 18(trT) 2 ) + Y? BY t (-6tvT + 12trB + 4tr£) 
~f 2 (-\ + k) + 7(1 + |)Ttr(T) + YlBY t {f b + ffc)] 

(9 - 3£;)T 2 + (27 - 9k) f trT + (9 - 3/^5^1 



+ <7l 



16(Jfe - l)TtrT + f (T 2 + y/sy t ) 



r,,,) 247 '' ll,4 tvr-f|trS-ftr£?] 



+ ghtn-L - l -§k) + ghln-f + fk) 



+ <#T(-f - ffc - 18ni - 6n 5 ) + g 2 2 g 2 3 T(-88 + 16k) 



T(f - 272jfe - 24n 10 - 8n 5 ) - f trT - f trS 



+ fifKf - i^io - 4^5 - §f& + §nio + fn 5 ) 

+ 4 2 (f - Sfcnio - #*n 5 - ™ k + f|n 10 + §n 5 ) 

+ </S(^ + f fcnio + fkn 5 + '-fk + 4n 10 n 5 + 2 -fn 10 + 6n 2 10 + 2 -f n , + | n J) 
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+ f2fs(60 - 4kn w - 12k + 20n 10 ) 



(2.20) 



7 g) = (Jfe + l) [3tr(S 3 ) + tr(E 3 )] + 9tr(T 2 S) + 18trTtr(TB) 

+ 6(tr£ + 3trS) [tr(E 2 ) + 3tr(S 2 )] + (24 - 8% 2 [tr(TB) + 3tr(S 2 )] 
+ g%[18ti(TB) + (9k + 9)tr(S 2 ) + (3k + 3)tr(£ 2 )] 

+ g 2 [-ftr(TB) + |fctr(TB) + 3tr(S 2 ) + f ktr(B 2 ) + 9tr(£ 2 ) - § ktr(E 2 )} 

+ 9$[-§trT - trB(m + ^ * + |Z„ 10 + g n5 ) + t r£(^/c - §f- (2.21) 

" f^io " §ri 5 )} + g 2 g 2 (-±trB - § HrB - %trE + %HrE) 

+ trB[g 2 g 2 (- 2 -§ + fffe) + g 2 g 2 (-132 + 24k)] + gl[-f trT 

- trS(^ + f k + f n 10 + fn 5 ) - trT(f + f /c + f mo + f n 5 )] 

- gttrB(±f + Ik + 48n 10 + 16n 5 ) + H 
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7 J 3) = (2k + 6)B 3 + 6YjT 2 Y b - 2Y h ] BTY b + Y^TY b (12trT - 6trS - 2tr£) 
- 2Y b ] TBY b + B 2 (6trB + 2tr£) 

+ S[12tr(TS) - 18(trS) 2 - 12trBtr£ + 36tr(S 2 ) - 2(tr£) 2 + 12tr(£ 2 )] 

+ y 2 y 6 t Ty 6 (-f + ffc) + <7 2 B 2 (-± + \k) + g 2 B[(7 - k)trB + (k - 3)trE] 

+ glY h ] TY b (9 - 3k) + g\B 2 (9 - 3k) + g 2 B[(27 - 9k)tiB + (9 - 3k)tiE] 

+ fglY^TYf, + fgjB 2 + # 2 B[(16/c - 16)trS + 16tr£] 

+ g\ [b(-*§ -JL.k- fn 10 - ln 5 ) - fftrT - ±§ tr B - §tr£ 

+ <fo 2 2 B(-f + Ifc) + ^ 3 2 5(-f + fk) + g 2 g 2 B(-88 + 16k) 



+ g%B(-f - § k - 18nio - 6n 5 ) 



B(f - 2|2 fe _ 24nio _ 8n5 ) _ 80 trT _ 80 tlB 



+ 9l( 2 -W + f *»io + f fens + ^ fe + 4n 10 n 5 + f> n 10 + 6n? + ™ n 5 + |n 2 ) 
+ ~ K^n 1Q + ^n 5 + if) + n 10 (In 5 + ±§ + > 10 ) + f§n 5 + ^n 2 ) 

+ fifKf - ^^io - iokn 5 ~^k + ±n w + ^n 5 ) 
+ atgiCM - i^io - EE*** - + ±f mo + £n 5 ) 
+ - l^io - - ff jfe + f§n 10 + fn 5 ) 

+ flf|(60 - 4/cnio - 12/c + 20ni ) 

(2.22) 



7 jg = 54trTtr(T 2 ) + 3(Jfe + l)tr(T 3 ) + 18trStr(TS) + 9tr(BTB) + 6trEtr(TB) 
+ g 2 i« 5 i- IkMT 2 ) + (! + J fc)tr(TS)) + <7 2 2 [(9 + 9fc)tr(T 2 ) + 18tr(TB)] 
+ gj[(72 - 24k)tr(T 2 ) + (24 - 8fc)tr(T£)] 

- [(W + mk + §n 5 + 4fn 10 )trT + §trB + %trE] (2.23) 

+ fe 2 (S^ - S) + alaK W k -'-¥) + ghlm - i32)]trT 

+ ^ 4 (trT(- 2 f - f k - f n 10 - f n 5 ) - f trB - f tr£) 

- ^(^trT + ffctrT + 48ni trT + 16n 5 trT) + S 
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7 ( 3 ) = (6 + 2k)E 3 + (6trS + 2tr£)£ 2 

+ S(12tr(TS) - 18(trS) 2 - 12trBtr£ + 36tr(S 2 ) - 2(tr£) 2 + 12tr(£ 2 )) 

+ glE\l + p) + g\E\9 - 3k) + &E{{lf + J fc)trS + (§ + p)trE) 

+ g%EtrB(16k - 64) + g%E(tiB(27 - 9k) + txE(9 - 3fc)) 

+ ^(-^ - gfc - f n 10 - fn 5 ) + g\g\E(-21 + fk) 

+ gK-f trT - f trfl - f tr£) + g%E(-f - § fc - 18n 10 - 6n 5 ) 

+ 9i C-M ~ f^io " ^5 - fi* + f n 10 n 5 + f n 10 + §n 2 + gn 5 + ^n 2 ) 

+ gigKf - & kn w - % kn $ - it k + jo n w + %n 5 ) 

,,,4 2/264 72t,_ 16^ 2647,, 72„ , 16„ \ 

+ 0103l~5 25 feni ~ 25 fen 5 ~ "25-^ + — ™10 + 

(2.24) 

where k = 6£(3), and 

•3 — «6/345 , l/45_ , 15 „ , 105 \ , „ /9_, , 81 , 27_ \ , 27„ , 3-,2\ 
= — 02l~ + fcl-g-^lO + -§"^5 + —J + ^10 1 4^5 + ~2 + "8" n 10j + "2"^5 + g n 5J 

1 „6/ 1839 j,/ 69 „ , 7 „ , 597\ , „ / 9 „ , 753 , 27„ \ , 211^ ,3 „2\ 
+ fllW - ^200^10 + 200 n 5 + 500~) + ™10l20™5 + TM + lo 71 ^ + 100 n 5 + 40%) 

+ 0102 (^j ~~ 2M^ n 10 - M) kn 5 ~ ~m k ~ M n 10 + M 77-5 ) 

+ 0i03 2 (f " i^io - 4*n 5 - ffe + fn 10 + |n 5 ) 



+ 0i02(f - igfcrcio - ~E k ~ T6 n io + fon 5 ) 

,4„2/ 



+ 020s( 90 - 6/cn io - 18/c + 30ni ). 

(2.25) 

In terms of the anomalous dimensions, the Yukawa /3-functions are: 

(3 Yt = lQ Y t + Y t ( lt + 7 h 2 ) , P Yb =l Q Y b + Y b ( lb + -y Hl ), Py t = 1lY t + Y T ( 7r + lHl ), 

(2.26) 

and the /3-function for the Higgs /U-term is 

P» = v(lH 2 +l Hl ). (2.27) 

We will also require the anomalous dimensions of the constituents of the extra 5 and 10 
representations, which are easily obtained by setting T = B = E = 0, except retaining 
terms that contain T,B,E only inside traces; such terms occur for the first time at three 
loops. 
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From the above expressions for f3 g . and 7 we have calculated the three-loop soft (3- 
functions using Eq. ( |2.3[ ) and FORM. The resulting expressions are very unwieldy; as an 
example we give the one, two and three-loop results for /3 m 2 , in the approximation that we 

Qt 

retain only g 2 and the top quark Yukawa coupling At (in what follows we denote the third 
generation squarks as Q t , t c , o c , and the first or second generation squarks as Q u , it c , d c ): 

/?« = 2\ 2 (Z t + A 2 t )-8(±g 2 M 2 + lg 2 M 2 + lg 2 M 2 ) (2.28a) 

Qt 

/?$ = -20A^(E t + 2A 2 t ) + lQglM 2 (n 5 + 3n 10 - §) 

Qt 

+ fgl(2m 2 Qt + m 2 c + m 2 c + (n 10 + 2)(m 2 uC + 2m 2 Q J + (n 5 + 2)m%) (2.286) 
0$ = [(1280* + + 16n§ + + 3 -fk)(n 5 + 3n 10 ) 

Qt 

+ 96n w n 5 + 144r^ )M 3 2 + (^ - f(n 5 + Sn w ))(m 2 c + m\ c + 2m 2 Qt ) 
+ (2m 2 Qu + m^ c )(f - f n 5 + fn 10 - fn 5 n 10 - 16n 2 ) 
+ m%(^f + ^n 5 - 32n 10 - 16n 5 n l0 - fnjM 

_ [(288 + ^fk + 48(n 5 + 3n 10 ))M 2 - (192 + ^k + 32(n 5 + 3n 10 ))A t M 3 
+ ( m k + m + 8{n5 + 3n 10 ))(Et + A 2 )]\ 2 g 4 3 
+ (lf + 32k) [M 2 - 2A t M 3 + E* + 2 A 2 } \jg 2 

+ (6k + 90)(E t + 3A 2 )X 6 t , (2.28c) 

where E t = + m| + m| c . For this special case, and also with n$ = nio = 0, the 
three-loop result, Eq. ( [2.284 ), was given in Ref. |2"2]| , except that in the corresponding 
expressions in this reference the squark masses of different generations are not clearly 
distinguished (as they must be since the third generation evolves differently from the other 
two). Complete results for the three- loop /3-functions including all three gauge couplings 
and n g x n g Yukawa matrices are available at Ref. . 

In our analysis we do include "tadpole" contributions, corresponding to renormalisa- 
tion of the Fayet-Iliopoulos (FI) .D-term at one and two loops. These contributions are 
not expressible exactly in terms of /3 9i ,7; for a discussion, and three-loop results for the 
MSSM, see Ref. ||23j| . For universal boundary conditions, the FI term is very small at low 
energies if it is zero at gauge unification; including the three-loop (FI) effects would have 
a negligible effect on our results. 
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3. The Snowmass Benchmark Points 

In this section we examine the effect of the three-loop corrections on the standard 
running analysis, that is for 715 = mo = 0. We will focus on the standard treatment with 
universal boundary conditions at gauge unification, often termed CMSSM or MSUGRA. 
Thus we assume that at Mx we have universal soft scalar masses (mo), gaugino massesl 

(mi) and A-parameters (A), and work in the third-generation-only Yukawa coupling ap- 

2 

proximation. This is for ease of comparison with existing results rather than because we 
find the scenario particularly compelling. We will present results for the set of MSUGRA 
Snowmass Benchmark Points shown in Table 1: 



Point 


tan (3 


mi 

2 


m 


A 


sign^ 


SPSla 


10 


250GeV 


lOOGeV 


-lOOGeV 


+ 


SPS16 


30 


400GeV 


200GeV 





+ 


SPS2 


10 


300GeV 


1450GeV 





+ 


SPS3 


10 


400GeV 


90GeV 





+ 


SPS4 


50 


300GeV 


400GeV 





+ 


SPS5 


5 


300GeV 


150GeV 


-ITeV 


+ 


SPS6 


10 




see footnote 


150GeV 





+ 



Table 1: Input parameters for the SPS Benchmark Points 



Other input parameters are shown in Table 2: 



pole 

mf 


pole 

m b 


pole 

171t 


a 3 (M z ) 


a 2 (M z ) 


oixiMz) 


178GeV 


4.9GeV 


1.777GeV 


0.1172 


0.033823 


0.016943 



Table 2: Input parameters for the running analysis 



In Table 2 the input couplings CK1...3 correspond to the Standard Model MS results; 
we calculate the appropriate dimensionless coupling input values for the running analysis 
by an iterative procedure involving the sparticle spectrum. We define gauge unification 

3 except for the SPS6 point. The SPS6 point corresponds to non-unified gaugino masses, 
Mi = 480GeV, M 2 = M 3 = 300GeV. 
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to be the scale where «2 and a.\ meet; we speed up the determination of this by (at each 
iteration) adjusting the unification scale using the solution of the one-loop /3-functions for 
the gauge couplings from the previous value of the scale. We employ one-loop radiative 
corrections as detailed in Ref. |24|] I; thus we run up from Mz using the full supersymmetric 
/3-functions. For most particles we evaluate the pole mass at a renormalisation scale equal 
to the pole mass itself, and determine this value by iteration; the exception being the light 
CP-even Higgs, where we use a scale equal to the average squark mass. 

3.1. Benchmark point SPS la 

This point is a "typical" point in MSUGRA parameter space. In Table 3 we compare 
our results for a selection of sparticle masses (at 715 = nio = 0) with the spread of results 
taken from Ref. ||, denoted AKP (note our convention that the predominantly right- 
handed top squark is £2)- 



In the first line of Eq. 37 of Ref. [24|, the first term in the square bracket should read 



— + rn 1 2 ) Bo(i r nj 2 , mi 1 , 0): i.e. it should have a minus sign. The corresponding exact result in 
Eq. D49 is correct, however. 
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mass 


Hoop 


21oops 


31oops 


AKP 


9 


628 


613 


611 


604 - 612 


ti 


594 


590 


583 


577 - 588 




400 


399 


391 


396 - 401 




573 


565 


557 


565 - 569 


ur 


552 


548 


539 


547 - 549 


h 


520 


514 


507 


514 - 518 


~b 2 


551 


548 


540 


539 - 548 


d L 


579 


571 


563 


571 - 574 


da 


551 


548 


539 


546 - 548 




212 


207 


206 


208 - 211 




139 


135 


135 


134 - 136 


eL 


209 


202 


202 


204 - 207 


CR 


147 


144 


144 


143 - 146 


i>e 


192 


186 


185 


186 - 191 


V T 


191 


185 


184 


185 - 191 


Xi 


104 


97 


97 


95.6 - 97.4 


X2 


193 


180 


179 


181 - 182 


X3 


351 


369 


364 


362 - 371 


X4 


376 


388 


384 


381 - 390 


4- 


193 


179 


178 


180 - 182 


+ 

X 2 


376 


388 


384 


380 - 390 


h 


114 


114 


114 


112 - 115 


H 


392 


403 


399 


403 - 407 


A 


391 


403 


399 


400 - 406 


H± 


400 


412 


408 


410-415 



Table 3: Sparticle masses (in GeV) for the SPSla point 



We would expect our two-loop results to correspond most closely to AKP and we see 
that they are indeed broadly consistent, typically being within the range defined by the 
other programs or within a GeV of it. The effect of inclusion of three-loop running is 
never greater than 2%; note, however, that the shift caused by three-loop running effects 
is comparable for ul and larger for t 2 ,UR than that produced by two-loop running effects. 
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3.2. Benchmark point SPS lb 

This is another "typical" point but with a higher value of tan j3. Our results are given 
in Table 4. 



mass 


Hoop 


21oops 


31oops 


AKP 


9 


967 


946 


943 


933 - 943 


ti 


848 


841 


832 


836 - 839 


t 2 


657 


656 


646 


652 - 661 


UL 


891 


878 


868 


878 - 882 


ur 


854 


849 


837 


848 - 850 


h 


781 


773 


763 


773 - 778 


h 


831 


827 


816 


819 - 828 


d L 


895 


882 


872 


882 - 885 


dn 


851 


847 


835 


844 - 848 


fi 


353 


347 


346 


347 - 349 




208 


199 


200 


196 - 202 


eL 


348 


339 


338 


341 - 342 


eR 


258 


254 


254 


253 - 256 


V e 


338 


329 


328 


329 - 332 


V T 


328 


318 


318 


319 - 322 


Xi 


173 


162 


162 


159 - 163 


X2 


327 


305 


304 


308 - 308 


X3 


507 


532 


526 


521 - 534 


X4 


526 


546 


541 


534 - 546 


xt 


327 


305 


304 


307 - 308 


xt 


526 


547 


541 


535 - 547 


h 


118 


118 


118 


117-119 


H 


528 


544 


539 


540 - 544 


A 


529 


545 


540 


538 - 544 


H± 


535 


551 


547 


547- 551 



Table 4-' Sparticle masses (in GeV) for the SPSlb point 
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3. 3. Benchmark point SPS 2 



This is a "focus point region" point ||25|| , characterised by the large value of m®. Our 
results are given in Table 5. 



mass 


Hoop 


21oops 


31oops 


AKP 


9 


835 


816 


814 


778 - 805 


ti 


1322 


1292 


1287 


1291 - 1318 


h 


942 


921 


913 


913 - 942 


UL 


1597 


1562 


1558 


1566 - 1591 




1584 


1556 


1552 


1556 - 1581 


bi 


1303 


1273 


1268 


1280 - 1309 


h 


1571 


1544 


1540 


1527- 1568 


d L 


1600 


1564 


1560 


1567 - 1593 


oIr 


1584 


1556 


1553 


1555 - 1580 


fi 


1463 


1454 


1454 


1455 - 1460 


T2 


1444 


1440 


1441 


1439 - 1443 




1468 


1459 


1459 


1460 - 1465 


en 


1457 


1453 


1453 


1453 - 1455 




1465 


1456 


1456 


1457 - 1463 


V T 


1459 


1450 


1450 


1451 - 1457 


Xl 


132 


123 


123 


121 - 124 


X2 


257 


237 


237 


240 - 241 


X3 


562 


579 


582 


528 - 596 


X4 


574 


589 


592 


539 - 605 


xt 


257 


237 


237 


240 - 241 


xt 


574 


590 


592 


539 - 605 


h 


119 


119 


119 


117-117 


H 


1548 


1545 


1546 


1542 - 1555 


A 


1548 


1545 


1546 


1532 - 1555 


H± 


1550 


1547 


1548 


1544 - 1557 



Table 5: Sparticle masses (in GeV) for the SPS2 point 
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3.4- Benchmark point SPS 3 

This is a "co-annihilation region" point, its distinctive feature being a light stau not 
much heavier than the neutralino LSP. Our results are given in Table 6. 



mass 


Hoop 


2 loops 


31oops 


AKP 


9 


964 


943 


940 


930 - 940 




851 


845 


835 


836 - 843 


h 


645 


644 


634 


640 - 650 


UL 


872 


860 


849 


861 - 863 


ur 


835 


830 


818 


828 - 831 


bi 


794 


787 


776 


786 - 793 


b 2 


830 


826 


814 


816 - 825 


d L 


876 


864 


853 


864 - 867 


dn 


831 


828 


816 


825 - 829 


fi 


300 


291 


290 


293 - 294 




180 


173 


173 


172 - 176 




299 


288 


288 


291 - 293 


£r 


186 


181 


181 


179 - 183 


T>e 


287 


277 


276 


277 - 281 


V T 


286 


276 


275 


276 - 280 


Xi 


172 


161 


161 


158 - 162 


X2 


325 


302 


301 


305 - 306 


X3 


512 


538 


531 


528 - 540 


X4 


533 


554 


548 


543 - 555 


xt 


324 


302 


301 


304 - 306 


xt 


533 


554 


548 


542 - 555 


h 


117 


118 


117 


116-118 


H 


579 


597 


591 


593 - 600 


A 


579 


597 


591 


589 - 600 


H± 


585 


603 


597 


598 - 605 



Table 6: Sparticle masses (in GeV) for the SPS3 point 
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3. 5. Benchmark point SPS4 



This is a point with large tan/3. Our results are given in Table 7. 



mass 


Hoop 


21oops 


31oops 


AKP 


9 


759 


743 


741 


729 - 738 


ti 


705 


700 


693 


693 - 697 


t 2 


544 


541 


533 


540 - 544 


UL 


777 


764 


757 


766 - 772 


ur 


755 


747 


739 


747 - 751 


h 


624 


619 


611 


614 - 619 


h 


693 


690 


683 


679 - 692 




782 


769 


761 


770 - 776 


oIr 


753 


746 


738 


746 - 749 


T\ 


423 


420 


420 


414-421 


T2 


272 


268 


268 


253 - 269 


eL 


455 


450 


449 


451 - 452 


eR 


419 


417 


418 


417-419 


l>e 


447 


441 


441 


442 - 445 


V T 


395 


390 


390 


387 - 393 


Xi 


128 


120 


120 


119 - 121 


X2 


242 


226 


225 


228 - 228 


X3 


400 


419 


415 


406 - 420 


X4 


420 


435 


431 


422 - 436 


xf 


242 


226 


225 


227- 228 


xt 


421 


436 


432 


422 - 436 


h 


116 


116 


116 


114 - 116 


H 


370 


386 


385 


355 - 367 


A 


371 


388 


387 


355 - 367 


#± 


381 


397 


396 


366 - 379 



Table 7: Sparticle masses (in GeV) for the SPS4 point 
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3. 6. Benchmark point SPS 5 



mass 


Hoop 


21oops 


Sloops 


AKP 


y 


743 

( j: u 


729 


727 


719 - 729 




653 

U U U 


654 

U U j: 


646 

U^tU 


629 — 651 

XJ^JiJ U U -L 


to 


265 

U U 


278 


263 


258 — 280 

Zj U O Zj OU 


It T 


684 


677 

U ( ( 


668 

UUO 


676 _ 685 

U I U UOU 




658 


656 


646 


655 — 660 


hi 


563 

U U U 


563 

U U U 


554 

U U j: 


U Ut: UU ( 


bo 


654 

U U j: 


653 

U U U 


643 


UUU uuu 


dr 


688 


681 


673 


681 _ 689 


u ti 


656 

U U U 


655 

U U U 


645 

V/ 1 U 


65S _ 658 

UUU uuo 


' 1 


264 


259 


258 


259 — 262 




186 


182 


183 

-LOO 


182 — 184 


P T 


263 


257 


257 


258 — 261 


P D 


195 

-L ' / U 


192 


193 


192 — 194 




251 


245 


245 


246 - 249 


it 


249 


243 


243 


244 — 247 


Al 


128 


120 


120 


119 — 120 




247 


229 


228 


230 — 236 

ZvUU A(UU 


Vo 

A3 


608 


626 


621 


626 — 631 

UajU UU-L 


v ,1 


621 


637 

U U i 


632 

U U £a 


6S7 _ 641 

UU ( Ut: J. 


Al 


247 


229 


228 


230 - 236 


X2 


620 


637 


632 


636 - 641 




117 


118 


118 


116-122 


H 


667 


682 


676 


681 - 694 


A 


667 


682 


677 


682 - 690 


H± 


672 


687 


681 


687 - 698 



Table 8: Sparticle masses (in GeV) for the SPS5 point with m t = 178GeV 

This point differs from the previous ones in having a large value of the A-parameter. 
The contributions of A to the off-diagonal term in the stop mass matrix have the same 
sign, and the magnitude of A is large, resulting in a light stop. For this point we have 
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calculated both using in Table 8 m t = 178GeV (as for the previous tables) and for com- 
parison in Table 9 with mt = 174.3GeV. This illustrates the sensitivity to the input mt, 
with the light stop changing over 20GeV due to this small change in mt. 



mass 


Hoop 


21oops 


31oops 


AKP 


9 


743 


729 


727 


718 - 728 


ti 


652 


653 


645 


628 - 649 


t 2 


243 


257 


240 


232 - 258 


u L 


684 


677 


668 


676 - 684 


ur 


658 


656 


646 


653 - 660 


h 


561 


560 


551 


551 - 564 


h 


654 


653 


643 


629 - 655 


dh 


689 


681 


673 


680 - 689 


d R 


656 


655 


645 


651 - 658 


T\ 


264 


259 


258 


258 - 262 


T2 


186 


182 


182 


182 - 184 


e-L 


263 


257 


257 


258 - 260 


eR 


195 


192 


192 


192 - 194 




251 


245 


245 


246 - 249 


V T 


249 


243 


243 


244 - 246 


Xi 


128 


120 


120 


119 - 121 


X2 


247 


229 


228 


230 - 236 


X3 


615 


632 


628 


632 - 637 


X4 


628 


644 


639 


643 - 646 


xt 


247 


229 


228 


230 - 236 


xt 


627 


643 


639 


643 - 646 


h 


115 


115 


115 


112 - 119 


H 


674 


688 


683 


687 - 693 


A 


674 


688 


683 


689 - 693 


#± 


679 


692 


687 


694 - 702 



Table 9: Sparticle masses (in GeV) for the SPS5 point with m t = 174.3GeV 
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3.7. Benchmark point SPS6 

This is a point with un-unified gaugino masses so we are unable to compare with 
Ref. |§. We instead use the paper by Ghodbane and Martyn (GM), Ref. PI , which also 
compares the results for various programs (Isajet, Susygen and Pythia). The results for 
these three programs are reasonably consistent with each other; this is due to some extent, 
however, to the fact that the Isajet gauge unification outputs are used as inputs for the 
other two programs; in our table we show only the Isajet predictions. Agreement with our 
results is less impressive; however we should notice that Ref. uses an earlier version 
of Isajet (7.58) than Ref. 0. Thus if we return to SPSla and compare the Isajet 7.58 
prediction for the gluino mass (595GeV) with the Isajet 7.69 one of 612GeV obtained from 
Ref. ||, we can anticipate that for SPS6 the more recent Isajet would give results more 
consistent with our (two-loop) ones, making the reasonable assumption that the newer 
version will give, for example, a higher gluino mass prediction for SPS6 as well. Our 
results for SPS6 are given in Table 10. 
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mass 


Hoop 


2 loops 


31oops 


GM 


9 


744 


726 


724 


708 


ti 


686 


681 


673 


661 


h 


498 


496 


488 


476 




684 


674 


665 


639 


ur 


665 


659 


650 


628 




620 


613 


605 


589 


h 


657 


652 


643 


624 


d L 


689 


679 


670 


644 


da 


658 


653 


644 


622 


Tl 


278 


271 


271 


270 


T2 


235 


229 


229 


228 


eL 


274 


266 


266 


265 


£r 


243 


238 


238 


237 


T>e 


261 


253 


253 


252 


V T 


260 


252 


252 


252 


Xi 


201 


190 


190 


189 


X2 


239 


222 


221 


218 


X3 


399 


419 


414 


399 


X4 


425 


439 


435 


420 


4- 

xt 


237 


220 


219 


215 


x 2 


423 


438 


434 


419 


h 


115 


115 


115 


115 


H 


469 


481 


477 


464 


A 


469 


481 


477 


463 


H± 


477 


489 


484 


470 



Table 10: Sparticle masses (in GeV) for the SPS6 point 

3.8. Discussion 

A clear feature of the results is that the corrections due to two and three-loop running 
can be quite large for squarks, but are typically smaller for weakly-interacting particles. In 
particular the light CP-even Higgs mass is very stable. The large three-loop a 3 corrections 
stem mainly from the M| contributions to the three-loop m 2 /3-functions; note that for 
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the only MSUGRA point such that uiq > mi, i.e. SPS2, the three-loop correction to the 

2 

squark masses is smaller than the two-loop one. 

Generally speaking we would anticipate that for regions of parameter space where the 
three-loop corrections are comparable to or exceed the two-loop ones, the four-loop ones 
will be at least as large. This suggests that we are already at three loops approaching the 
asymptotic region for the /3-functions. So it appears that squark mass predictions with an 
accuracy greater than a few per cent will not be possible using perturbation theory. 

Overall our results agree reasonably well with those of existing programs ||. One 
place where we have a significant difference is for the H, A, results for SPS4. This is 
a large tan/3 point; however our results for the 6-squark and (i-squark masses (which one 
would expect to be sensitive to large tan/3) agree quite well, so for the moment we have 
no explanation for this discrepancy. 



4. The Semi-perturbative Region 

The addition of additional matter representations in complete SU5 multiplets does 
not affect gauge unification (and the unification scale) at one loop. Beyond one loop this 
is no longer the case, and increasing the amount of matter relevant to the running analysis 
requires the presumption of larger threshold corrections at the unification scale in order to 
restore gauge unification; one is thus forced to argue that the success of gauge unification 
in the MSSM is coincidental @. 



5 Historically gauge unification was implemented by using a^(Mz) as an input and computing 
sin 2 6w, although the latter was more accurately measured, because sin 2 9w varies very slowly 
with a-i{Mz)i and conversely (of course) 03 (Mz) varies rapidly as a function of sin 2 #w The 
current experimental results for 0:3 (Mz) already require us to suppose the existence of some high 
scale radiative corrections in the MSSM; but the fact remains that things get worse as we add 



more matter [ 13 j. 
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a 1 ,a 2 ,a 3 vs x for SPS1 a 
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Fig.l: Gauge coupling unification for n±o = 1.7. Solid, dashed, and 
dotted lines correspond to a\,ai,a^ respectively. 

In Fig. 1 we show the evolution of the gauge couplings ai = gf/^Ti) for nio = 1.7, 
using three-loop /3-functions for all couplings. (As remarked in Ref. [0, the mass scale 
of these additional multiplets being unknown it makes sense to parametrise their effects 
by taking 77,5,77-10 to be continuous variables.) The couplings are plotted against r = 
2^: ln((3/Mz); evidently we are still in the perturbative regime. The input parameters at 
Mz correspond to a typical supersymmetric mass spectrum; specifically, the Benchmark 
point SPSla. One sees clearly the need for large corrections to restore gauge unification. 

We gave a number of examples of the effect of additional matter on the sparticle 
spectrum predictions in a previous paper [lf|; here we contrast the effect on the first and 
third generation squark masses. Thus in Fig 2 we plot, for the SPS5 point, the ratio of 
the ul and gluino masses against nio for 775 = 0; as already noted in Ref. [^2[ , the mass 
increases with 7710. It is interesting that the effect of the three-loop correction to this ratio 
almost precisely cancels the two- loop correction, for all nio- We contrast this with Fig 3 
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where we show the behaviour of the light stop mass for the same SPS point; in this case 
the ratio decreases smoothly, and the three-loop correction only cancels the two-loop one 
at nio = 0. For the SPS5 point the electroweak vacuum fails around nio = 0.48. (The 



change in this value and in Fig. 3 from our previous paper [15] is due to the change in the 



input top pole mass, and to an improved treatment of the Higgs potential minimisation.) 

In Fig. 4 we plot the light CP-even Higgs mass for SPSla as a function of nio (for 
715 = 0). We see that it is fairly stable both with respect to loop corrections and the 
addition of extra matter. In the case of SPSla the electroweak vacuum fails at around 



nio = 1-8. 



1.01 



u L /gluino mass ratio vs n 1Q for SPS5 



0.98 - 




0.92 



0.91 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Fig. 2: Plot of the u^jgluino mass ratio against nio for SPSS. Solid, 
dashed and dotted lines correspond to one, two and three-loop running 



respectively. 
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0.4 



Light stop/gluino mass ratio vs n for SPS5 




0.3 - 



! 0.25 - 



0.2 - 



0.1 5 ' 1 1 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

n io 

Fig. 3: Plot of the light stop/gluino mass ratio against nio for SPS5. 
Solid, dashed and dotted lines correspond to one, two and three-loop 
running respectively. 

Light Higgs mass vs n for SPS1a 

117 | 1 1 1 1 1 1 1 1 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

"to 

Fig. 4: Plot of the light CP-even Higgs mass against nio for SPSla. 
Solid, dashed and dotted lines correspond to one, two and three-loop 
running respectively. 
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5. Conclusions 



We have extended typical detailed running coupling analyses for the MSUGRA MSSM 
SPS benchmark points to incorporate three-loop /^-function corrections for the running 
masses and couplings. We compare our results to those obtained by existing programs 
using two-loop running. The spread in the results from these programs is probably due 
to a mixture of program errors and genuine theoretical uncertainties such as the choice of 
scale appropriate for the evaluation of the pole mass. Presumably over time the results 
used by these programs will converge; we would argue that a more reliable estimate of 
the ultimate theoretical error in these spectrum calculations is currently provided by the 
difference between our two and three-loop calculations, as opposed to the spread in the 
various available two-loop results. 

Generally speaking the effect of the three-loop running corrections is small for weakly- 
interacting particles but larger for the squark masses. For the light stop mass at the SPS 5 
point, we see an 8% effect, but more typically the effect is between 1% and 2%. This 
appears to us to represent a fundamental limit on the theoretical precision of squark mass 
theoretical predictions. 

Finally we show how additional matter in SU5 multiplets can affect the sparticle spec- 
trum; more dramatically as the "semi-perturbative unification" regime [ff2| is approached. 
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